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Abstract In the single degenerate (SD) scenario for type la supernovae (SNe la), a 
mass-accreting white dwarf is expected to experience a supersoft X-ray source (SSS) 
phase. However, some recent observations showed that the expected number of the mass- 
accreting WD is much lower than that predicted from theory, whatever in spiral or ellip- 
tical galaxies. In this paper, we did a binary population synthesis study on the relative 
duration of the SSS phase to their whole mass-increasing phase of WDs leading to SNe 
la. We found that for about 40% progenitor systems, the relative duration is shorter than 
2% and the evolution of the mean relative duration shows that it is always smaller than 
5%, whatever for young or old SNe la. In addition, before SNe la explosion, more than 
55% progenitor systems are experiencing a dwarf novae phase, and only no more than 
10% is staying SSS phase. These results are consistent with the recent observations, and 
imply that both in early- and late-type galaxies, only a small fraction of mass-accreting 
WD resulting SNe la contribute to the supersoft X-ray flux. So, although our results are 
not directly related to the X-ray output of SN la progenitor, the low supersoft X-ray lumi- 
nosity observed in early type galaxies may have no ability to exclude the validity of SD 
model. On the contrary, it is evidence to support the SD scenario. 
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1 INTRODUCTION 

Although Type la supernovae (SNe l a) are very impo r tant in astrophysics, i.e . as standard candles 



to measure cosmological para metersdRiess et al. 19981: iPerlmutter et aL 19991). t he exact nature of 



their progenitors is still unclear dHillebrandt & Niemeverl 20001: iLeibundgutl 20001: iParthasarathv et all 
120071) . There is a consensus that SNe la result fro m the thermonuclear explosion of a carbon- 
oxygen white dwarf (CO WD) in a binary system dHovle & Fowled Il960l) . According to the na- 
ture of the companions of the mass accreting WDs, two basic scenarios for the progenitors of 
SN la have been discussed over the last three decades. O ne is the single degenerate (SD) model 
dWhelan & Ibenlll973b iNomoto. Thielemann & Yokoilll984l) . i.e. the companion is a main-sequence 
or a slightly evolved star (WD+MS), a red giant star (WD+RG) or a helium star (WD + He 
star) dLi & van den Heuvellll997l lHachisu et alJI 19991 [Langer et al.ll2000l: lHan & Podsiadlowskill2004 : 
Chen & Lill2007t iMeng. Chen & Had 120091: iLii et all 120091: IWang et alj|2009t iMeng & Yandl2010c : 
Wang. Li & Hadl2010h. The other is the double degenerate (DD) model, i.e. the companion is another 



CO WD dlben & Tutukovlll984IWebbin9l984l) . 

It is believed that the X-ray signature of these two possible scenarios is very different, i.e. no strong 
X-ray is expected from the DD scenario until shortly before SNe la explosion, while the mass-accreting 
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WDs may show the signatures of supersoft X-ray sources (SSSs) long befor e it explodes as SNe la. This 
provides a potential method to distinguish the SD and DD model. Recently, iGilfanov & Bogdanld2010b 
(hereafter GB10) obtained the 0.3-0.7 keV soft X-ray luminosity and the K-band luminosity of several 
early type galaxies. By comparing the detected X-ray luminosity with those predicted from the SNe la 
rate estimated from the K-band luminosity, they found that the observed X-ray flux from the early type 
galaxies is a factor of ~ 30 - 50 less than predicted from the SD accretion scenario, and then they con- 
cluded that no more than five percent of SNe la in early t ype galaxies can be p roduced by mass-accreting 



white dwarfs of the SD scenario. However, as noticed by Hachisu et al.1 (1201 Ob . the results of GB 1 based 



on an assumption that all the accreting WDs are in the SSS phase lasting for about two million years 
before SNe la explosio n. This assumption is rather simple and may overestimate the duration of SSS 
phase by a factor of 10. lHachisu et al.l (1201 Ob recalculated their symbiotic model and found that the re- 
sults in GB10 str ongly supports the SD scenario as a main contributor of SNe la in early type galaxies. 
Similar to GBlO. lDi Stefanol(l2010l) also noticed that the predicted SSS number from SD model in spi- 
ral galaxies is also much larger than that from observations, which may imply that the d uration of SSS 
phase from SD scenario is much shorter than the accreting time of WD, as presented in lHachisu et al.l 
( 12010b . Actually, many works have referred to this problem, i.e the duration of SSS phase for an ac- 



creting WD is very short (s ee lYungelson et al 



1996i lHan & Podsiadlowskll2004l iMeng. Chen & Hani 



l2009HMeng & Yandl2010ah . In lMeng & Yangl fcOlOal). thev calculated a dense model grid including 



WD + MS and WD + RG systems, and the calculations obtained a very short duration of SSS phase for 
an accreting WD , i.e. th e typical value is 1.5 — 4.5 x 10 5 years whic h is consistent with the calculation 
of lHachisu et aL ( 2010b (see figures 1 and 2 in lMeng & Yandl2010ab . In this paper, we want to use the 



calculations in lMeng & Yang ( 2010ab t o do a binary population synthe sis study to check whether or not 
the SD model in lMeng & Yangl d2010ah may explain the discoveries of lDi Stefanold2010b and GB10. 



In section |2j we simply describe our method, and present the calculation results in section [3] In 
section|U we show discussions and our main conclusions. 

2 METHOD 

Recently, IMeng & Yangl d2010ab constructed a comprehensive single deg enerate progenitor m odel for 
SNe la. In the model, the mass-stripping effect of optically thick w i nd (lHachisu et ail 1 1996b and the 
ef fect of a thermally in stable disk were included (lHachisu et alj|2008t IXu & Lill2009b . The prescription 



of lHachisu et alj d 1999b on WDs accreting hydrogen-rich material from their companions was applied to 
calculate the WD mass growth. The optically thick wind and the material stripped-off by the wind were 
assumed to take away the specific angular momentum of the WD and its companion, respectively. In 
lMeng& Yangl d2010ab . both WD + MS channel and WD + RG channel are considered, i.e. Roche lobe 
overflow (RLOF) begins at MS or RG stage. The Galactic birth rate of SNe la derived from that model is 
comparable with that from observations. In a ddition, that model ma y even explain some supernovae with 
low hydrogen mass in their explosion ejecta (IMeng & Yang|2010bb . In their calculations, the progenitor 
system may experience a wind (Case wind), a calm hydrogen-burning (Case Calm), a nova explosion 
(Case Nova) or a dwarf nova phase (Case DNova). If a system is undergoing a calm hydrogen-burning 



phase, it may show the properties of SSSs. IMeng & Yang (1201 Oab calculated more than 1600 WD close 



binary evolutions and showed the initial parameter space leading to SNe la in an orbital period - sec- 
ondary mass (log Pj, Mg) plane. Their calculation is comprehensive, and the evolution details of all the 
models are obtained but not sorted for publishing. In this paper, we extract the evolutionary properties 
of the models from the data f iles of the calculations and incorporate them into the BPS code developed 
by iHurlev et al.l (120001 120021) to obtain the relative duration of SSS phase to the whole mass-increasing 

time of WD leading to SNe la. 

We simply describe the assumptions for mass accretion in lMeng & Yangld2010ab . In an initial binary 
system, the companion fills its Roche lobe at MS or during HG or RG, and mass transfer occurs. They 
assume that the transferred materials form a disk around the CO WD. If the mass-transfer rate is lower 
than a critical mass-transfer rate, the disk is thermal unstable and a dwarf nova is expected, otherwise the 
WD accretes the transferred material smoothly at a rate M a = |M tr |, where the critical mass-transfer 
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rate is a function of the orbital period of binary system: 



M C) th = 4.3 x l - 9 (^) 1 - 7 Af Q yr- 1 , (1) 



where P or b is the orbital period unit in hour (lOsakill 19961 Ivan Paradusll 19961). If the accretion r ate of 



CO WD, M a , is larger than a critical value, an optically thick wind occur (IHachisu et al 1 11996b . where 
the critical accretion rate is 

M C)H = 5.3 x 1(T 7 (L7 ~ X) (Mwd - O.4)M yr-\ (2) 
A 

where X is hydrogen mass fraction and Mwd is the mass of the accreting WD (mass is in Mq and 
mass-accretion rate is in M Q yr" - 1 . IHachisu et~aHll999h . If M a is smaller than M c h, but higher than 
|M Cj h, the hydrogen-shell burning is steady, i.e. the system is experiencing a SSS phase. When M a is 
lower than |Mc,H but higher than |M C ,H> a very weak shell flash is triggered and then a recurrent nova 
is expected. When M a is lower than gM C) H> the shell flash is so strong that no material is accumulated 
on to the surface of the CO WEQ. 

We extract the e volutionary time of Case calm from all the binary system calculated in 
lMeng& Yangl(l2010ah and found th at the typical value fo r SSS phase is 1.5 — 4.5 x 10 5 years which 
is consistent with the calculation of IHachisu et al.l d2010l) . Since IHachisu et all d2010t) have discussed 
it, we will focus on the relative duration of SSS phase to the whole mass-increasing time of CO WD 
until SNe la explosion, isssAsN, i- e - the fraction of the time spent in the SSS phase. Please keep in 
mind that as observational X-ray output of a nuclear-burning WD nonlinearly increases with its mass, 
i.e. it depends on the evolutional details of WD binary systems, the results obtained in this paper are not 
directly related to the X-ray output of SN la progenitor. However, the results here may partly reveal the 
X-ray nature of SN la progenitor at least. 

To obtain the distribution of tsssAsN, w e have performed a series detailed Monte Carlo simula- 
tions via Hurley's rapid binary evolution code dHurlev et al.ll2000ll2002l) . In the simulations, if a binary 
system evolves to a WD + MS or WD + RG stage, and the system is located in the (log P\ M\) plane 
for SNe la at the onset of RLOF, we assume that a SN la is produced, tgss / i$N is then obtained by 
interpolation in the th ree-dimensional grid (M ^ m , M\, log P 1 ) of the more than 1600 close WD binary 
ysytems calculated in lMeng & Yand d2010al) . In t he simulations, we foll ow the evolution of 10 7 sampl e 
binaries. The evolutional channel is described in iMeng & Yand" feOlOal) . As in lMeng & Yangl d2010al) . 
we adopted the following input for the simulations. (1 ) A single starburst or a constant star formation 



rate are assumed. (2) The initial mass function (IMF) of Miller & Scald (1 1979b is adopted. (3) The mass- 
ratio distribution is taken to be constant. (4) The distribution of separations is taken to be constant in 
log a for wide binaries, where a is the orbital separation. (5) A circular orbit is assumed for all bina- 
ries. (6)The common envelope (CE) ejection efficiency qp E, which denotes the fr action of the released 
orbital energy used to eject the CE, is set to 1.0 or 3.0. (See Mens & Yang 2010a for details). 



3 RESULTS 

3.1 the relative duration of SSS phase 

In figure [Tj we show the current-epoch-snapshot distribution of tsssAsN f° r different acE, where a 
constant star formation rate is assumed. We can see from the figure that the distribution is similar to 
each other for different the acs- A remarkable feature in the figure is that it is about 40% of mass- 
accreting CO WDs for SNe la that spend less than 2% of their mass-increasing time on SSS phase. 
Actually, for most of these WDs, the relative duration of SSS phase is only about 10~ 3 — 10~ 4 . No WD 
has a relative duration of SSS phase more than 30%, which means that if we assume that the WDs show 
the signatures of SSSs during their whole mass-increasing time until SNe la, the number of SSSs is at 



Actually, the mass of WD is reduced when a WD undergos a nova explosion. 
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Fig. 1 Distribution of the relative duration of SSS phase to the whole mass-increasing time 
of CO WDs until SNe la explosion, isssAsN, for constant star formation rate with different 
(*CE (solid nne: "ce = 10; dashed line: «ce = 3.0). For acE = 3.0, the line is up-moved 
by 0.1. 



leat overestimated by a factor of 3. In the distributions, there is a smaller bump between 0.16 and 0.2, 
which correspond to the syste ms with an initial WD mass of 0.6 ~ 0.8M Q , or a primordial primary 
of 2.5 ~ 3.5M Q (see figure l2l iMerig et al.ll2008l) . In figure [2] we give the distributions of the relative 
duration of SSS phase and the initial WD mass for qce = 1.0 (The distributions of qce = 3.0 is similar 
to that of acE = 1-0). From the figure, we can see that whatever the initial WD mass is, the low relative 
duration of SSS phase (less than 2%) is dominant. 

To explore the contribution of the accreting WD on the soft X-ray flux of galaxies, a mean relative 
duration of SSS phase is necessary. Since the relative duration of SSS phase is covering two orders of 
magnitude, we use a geometric mean to denote their mean value. Figure [3] show the evolution of the 
mean relative duration of SSS phase for a single star burst. It is amazing that the mean relative duration 
is always smaller than 5% from ~ 10 8 — 10 10 yrs, which may imply that even in spiral galaxies, only a 
small fraction of the mass-accreting WDs leading to SNe la contribute to supersoft X-ray flux. 

3.2 class 

In the calculations of iMeng & Yand d2010al) . the SNe la may explode at wind, SSS, recurrent nova 
and dwarf nova phase. In figure |4] we show the relative importance of wind, SSS, nova and dwarf 
nova at SNe la explosion. We can see from the figure that more than 55% progenitor systems show 
the properties of dwarf nova before SNe la, and about 30% are recurrent novae. Only no more than 
10% progenitor systems are experiencing the SSS phase. There are about 2% -10% (depending on cice) 
SNe la exploding at wind phase. These SNe la may show the properties of 2006X-like SNe la (see 
IMeng. Yang & Geng|l2010l) . Then, even before SNe la explosion, SSS is still not a dominant state. 
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Fig. 2 Distribution of the relative duration of SSS phase and initial WD mass with qce = 1.0. 



Figure [5] gives the evolution of the fraction of SNe la exploding at wind, SSS, recurrent nova and 
dwarf nova phase for a single starburst. Here, we only present the case of acE = 1.0 since the case of 
«ce = 3.0 is similar. We can see from the figure, when age is longer than 2 Gyr, there is no SNe la 
exploding at SSS phase and almost all SNe la exploding at dwarf nova phase. C onsidering that supersoft 
X-ray radiating from dwarf nova is rare (see discussions in lHachisu et al l2O10b . the discovery by GB 1 0, 
i.e. the supersoft X-ray flux is very low in early type galaxies, is a natural result. 



4 DISCUSSION AND CONCLUSIONS 

In this paper, we give the distribution of the relative duration of SSS phase and the evolution of the dis- 
tribution. Our results show that the duration of the SSS phase for accreting WDs is very short compared 
with their mass-increasing time. This result imply that if we assume that WDs radiate supersoft X-ray 
during their whole mass-increasing phase, the number of SSS should be overestimated significantly. The 
statement in GB10 that no more than 5% of SNe la in early type galaxies can be produced by mass- 
accreting white dwarfs of the SD scenario is just based on this simple assumption. Our BPS study shows 
that the mean relative duration of SSS phase is always smaller than 5% for SNe la with bo th short- and 
long-delay time, which is consistent with the discoveries of GB10 and iDi Stefanol d2010h . Our results 
mean that in spiral or elliptical galax ies, the fraction of mass-accreting WD contributing t o supersoft 
X-ray is very small, as suspected by |5i Stefanol d2010b . So, the discoveries of GB10 and IDi Stefanol 
(|2010|) could be evidence to support the SD scenario as SNe la progenitor model, and then we can not 
use the low supersoft X-ray luminosity of galaxies or the low number of mass-accreting WD showing 
supersoft X-ray signal in elliptical galaxies to exclude the validity of SD model. 

In addition, please keep in mind that there is another uncertainty in GB10 that they assume all the 
supersoft X-ray come from the mass accreting WDs which will explode as SNe la finally. Actually, 
the SSS at t = 10 10 yr are most likely not the progenitor of SNe la, whatever in spiral or elliptical 
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Fig. 3 Evolution of the mean relative duration of SSS phase to the whole mass-increasing 
time of CO WD until SNe la explosion for a single starburst with different «ce (solid line: 
otQE = 1-0; dashed line: acE = 3.0). 



galaxies dYungelsonll20 1 0t) . For example, some progenitor systems may experience a SSS phase, but 
nova ex plosion prevents them from becoming SNe la finally (see the figures 3 and 4 in iMeng & Yand 
d2010al) 1. Then, the contribution of mass-accreting WDs r esulting in SNe la to supersoft X-ray flux may 
be smaller than those shown in this paper, as discussed in lHachisuetal1d2010h . 

However, there are some uncertainties for our results in this paper. Firstly, as well known, a WD 
radius decreases as its mass increases. Correspondingly, at the same accretion rate M and luminosity , 
less massive WD would have lower temperature, and then more soft X-ray dvan den Heuvel|[l992 ). 
Calculation show s that WDs less massive than about 1 Mq do not contribute to the statistics of SSS 
dDi Stefandlioioh . The soft X-ray emission from galaxies is dominated by the emission from most 
massiv e WDs (maybe close to Ch andrasekhar mass limit because of the special X-ray band, i.e. 0.3-0.7 
keV in lGilfanov & BogdarfeOlOl) . For the reason above, the total luminosity of SN la progenitors not 
only simply scales with the duration of the super-soft phase but also depends on the details of their 
evolutional tracks. Si nce the most of SSSs happen when WD masses are around 1.2 — 1.3M© (see 
IMeng & Yan3l2010ah . it should be taken as an upper limit of the real cases for the number of SSSs 
deduced from the results in this paper. Then, this uncertainty can not change our basic conclusion, i.e. 
GB10 overestimates the number of SSSs. 

Secondary, a part of accreting ONeMg WDs may also contribute to the soft X-ray and their contri- 
bution was not considered here, which may make us to underestimate the number of SSSs. Fortunately, 
the number of ONeMg WDs must be much smaller than that of the CO WDs considered here because 
of initial mass function. So, the basic conclusions in this paper are still hold. 

Thirdly, in this pap er, we consider a new possible mechanisms to SNe la, i.e. dwarf nova, which is 
from a disk instability (Osakil ll996t Ivan Paradusll 19961) . For a low amount of accreted materials during 
an outburst, the accreted matter may simply accumulate onto the WD but remains unburnt. After a few 
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Fig. 4 The relative importance of wind, SSS, nova and dwarf nova at SNe la explosion for a 
constant star formation rate with different «ce (solid line: «ce = 1-0; dashed line: «ce = 
3.0). 



tens of disk instabilities, the accreted materials finally reaches a critical mass which may result in a shell 
flash. This could be nothing else but a classical (or recu rrent) nova (see also lHachisu et al.ll2010h . This 
scenario was uphold by observations (Wils et a l.luOlOl). Ho wever, its SSS duration is too short to be 
compared with the steady SSS duration (lHachisu et alj|2010l) . Then, our result, i.e. the dwarf nova is the 
dominant progenitors for old SNe la, may imply that the total luminosity of SSSs in early type galaxies 
would be much smaller than that predicted b y GB10, which is consisten t with the observation of GB10. 
In addition, I have to remind the reader that iGilfanov & Bogdanl d201 lb recently argued that WD with 
unstable nuclear burning, i.e. classical/recurrent nova, may not be the main class of SN la progenitor 
in early type galaxies because they would be producing too many recurrent nova outbursts inconsistent 
with classical nova statistic in nearby galaxies, i.e. the observational rate of classical nova with decay 
time shorter than 20 days is overestimated by a factor of 45-75. This inconsistency seems to be another 
evidence against the validity of SD scenario in early type galaxies. Actually, in this paper, the dwarf 
no va is the dominant one for the progenitor of SNe la in early type galaxies based on the calculation 
of IMeng & Yand d2010al) . In IMeng & Yand d2010ah . the duty cycle for dwarf nova is set to be 0.01 
(typical range is from 0.1 to a few 10~ 3 ), which means that if the d warf nova is a possible chan nel for 
SNe la in early type galaxies, the estimated rate of classical nova in IGilfanov & Bogdanl (1201 ll) could 
be overestimated by a factor of 100 (typical range is from 10 to several 10 2 ). Considering that there are 
sti ll some recurrent nova contr ibuting to SNe la in early type galaxies, the overestimated factor of 45-75 
bv IGilfanov & Bogdanl d201 ll) could be a natural results and should also not be a key against the SD 
channel in early type galaxies. 

Finally, some mass-accreting WDs may not explodes as SNe la and the contribution of these WDs 
to soft X-ray luminosity was not considered i n this paper. However , parameter space for these WDs is 
much smaller than those exploding as SNe la dMeng & Yandl20~10al) . Furthermore, the relative duration 
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Fig. 5 Evolution of the fraction of wind, SSS, recurrent nova and dwarf nova at SNe la 
explosion for a single starburst with qce = 1-0. 



of SSS phase for these WDs is much smaller than those exploding as SNe la. So, the missing WDs may 
not significantly affect our discussions in this paper. 
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